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1. INTRODUCTION 


Tue theory of the specific heats of crystals developed in a recent paper by 
the present author is based on the result that the structure of a crystal con- 
taining p atoms in each unit cell has (24p — 3) normal modes of vibration 
characterised by sharply defined monochromatic frequencies. These may be 
described as the modes of internal vibration of the groups of 8p atoms each 
included in volume elements of the crystal whose dimensions are twice as 
large in each direction as those of the unit cell. The three degrees of freedom 
of atomic movement left out in this description represent the three translations 
of each group of 8p atoms. The vibration spectrum resulting from these 
translations may be designated as “the residual spectrum ”’ and its nature 
was discussed in the paper already cited. The contribution to the thermal 
energy of the crystal arising from the translations of each group of 8p atoms 
is represented by an integral in which Einstein’s expression for the average 
energy of an oscillator of frequency v appears multiplied by a quantity 
9. v*dv/v> and the product is integrated between the limits 0 and 1. The 
upper limit is the upper limit of frequency in the residual spectrum and is 
also the lowest of the (24p — 3) frequencies characteristic of the structure. 


In the case of the metallic elements now under consideration, p = 1 
and the residual spectrum thus represents only 3 out of every 24 degrees 
of atomic freedom of movement. It nevertheless plays a major role in 
determining the magnitude of the specific heat of the metals at very 
low temperatures. Hence, a comparison of the experimental data at such 
low temperatures with the calculated values is of importance for establishing 
the correctness of the present approach to the theory of their specific heats. 
It is the object of the present paper to present the results of such a com- 
parison. 

As the determinations by Giaque and his collaborators which have so 
far been the basis of our discussion stop at 15° K., we have necessarily to 
make use of the results reported by other investigators for still lower tempera- 
tures. The circumstance that the specific heats at such temperatures are 
very low and increase rapidly with rise of temperature sets a limit to the 
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possible accuracy of their experimental determinations. The specific heat 
also varies enormously from metal to metal and hence even a small per- 
centage of impurity would affect the determined values, raising them or 
lowering them, as the case may be. These factors may be responsible for 
the divergences noticed between the results reported by Giaque and col- 
laborators and those obtained by other investigators in the temperature 
ranges common to them. These divergences are particularly large in the 
case of aluminium which has a very low specific heat in the range 0 to 
15° K. but are less serious in the case of the other three metals, copper, 
silver and lead. 


2. THe ELECTRONIC SpEciIFIC HEATS 


The observed variation of the specific heats of metals with temperature 
in the vicinity of the absolute zero differs in its features from that observed 
in the case of insulators, e.g., diamond or the alkali halides. The curves 
do not touch the axis of the abscisse but meet it at the origin making a definite 
angle with it. The specific heats in this region can be represented by an 
expression of the form AT + BT*. These facts are explained on the basis 
that the first term represents the contribution arising from the thermal agita- 
tion of the electrons and the second from the thermal agitation of the atoms 
themselves. This interpretation finds support in the fact that in the case of 
metals exhibiting superconductivity, e.g., aluminium and lead, the co- 
efficient A exhibits a distinctly larger value in the range of superconductivity, 
while the coefficient B remains unaltered. 


TABLE I 
Specific Heats at Low Temperatures 
Cc, = AT + BT® 
| | 
Range of | 
Substance measure- Ax 104 Bx 10° | Author 
ment | | 
| | | | 
Aluminium ..;  0-20°K. 3:48 | 0°63 | Kok & Keesom | 
0- 4° K. 3-42 | 0-68 | Howling et al. 
Copper ii 0-20° K. | 1-78 | 1-24 | Kok & Keesom 
0- 5° K. | 1-64 1-14 Corak et al. 
Silver ai 0- 5° K. 1-46 4-06 Corak et al. 
Lead “a 0- 5° K. 7-48 52-04 Horowitz et al. 
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Table I shows the results for the four metals now under consideration 
reported by various authors. It will be seen that the coefficient A is largest 
for lead; aluminium stands next and is followed by copper and silver, the 
coefficients for the latter not being very different. The value quoted in the 
case of lead is for the state of the metal in which the superconductivity has 
been eliminated by a strong magnetic field. 


It follows from what has been stated above that for a comparison of 
the experimental data with the results of the present theory, it is necessary 
to add to the specific heat a term representing the electronic contribution. 
This has been done, making use of the results of the authors quoted in the 
table. We shall discuss first the results thus obtained in the case of copper, 
silver and lead, since in the cases of these metals the determinations which 
we make use of are in reasonable agreement with those reported by Giaque 
and collaborators. The case of aluminium where this is not the case will 
also be discussed. 


3. COMPARISON OF THEORY AND EXPERIMENT 


Copper.—The specific heat of copper in the range between 0° and 20° K. 
has been reported by Kok and Keesom, while in the more limited range 
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from 0 to 5°K., measurements have been made by Corak, Garfunkel, 
Satterthwaite and Wexler. We present the comparison between theory and 
experiment in the form of two graphs, one in the range 0 to 20° K. and the 
other in the range 0 to 5°K. In the latter range the results of Kok and 
Keesom are systematically in excess of those reported by Corak et al., the 
difference being about 8% at 5°K. However, a satisfactory overall agree- 
ment between theory and experiment is exhibited by both the graphs. 


w 0°003 

w COPPER 

= oO OBSEVATIONS ANO KEESOM 

50-002F 

= 

a 

w 

=x 

= THEORETICAL CURVE INCLUDING 
s ELECTRONIC SPECIFIC HEAT 

i 


! 2 3 4 5 
ABSOLUTE TEMPERATURE 


Fic. 2. Specific Heats of Copper. 


Silver —Likewise in the case of silver we present the results in the form 
of two graphs, Figs. 3 and 4 respectively. Here again, a satisfactory agree- 
ment emerges between the theoretical results and the experimental data. 


Lead.—All the four characteristic frequencies besides the residual 
spectrum contribute to the observed specific heat at 20° K. and this con- 
tinues to be the case upto 10° K. At still lower temperatures the contributions 
of only the lowest of the four characteristic frequencies and the residual 
spectrum predominate. Below 5° K., only the residual spectrum is of import- 
ance. As in the previous cases we present a comparison between theory 
and experiment in the form of two graphs covering the ranges 0 to 20° K. 
and 0 to 5° K. as Figs. 5 and 6 respectively. In Fig. 5 the available determina- 
tions of various authors have all been plotted and it will be seen that between 
15° and 20° K. the determinations agree generally with each other and with 
the course of the theoretical curve. Between 5° and 15° K. the only observa- 


tions are those by Keesom and Van den Ende whose results appear to be 


rather erratic but nevertheless follow the trend of the theoretical curve fairly 
well. 


At the lowest temperatures the results of these same authors are notice- 
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ably in excess of those deduced from theory. But as will be seen from 
Fig. 6 the observations by Horowitz ef a/, are in excellent accord with the 
calculated values, “ 
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Aluminium.—Of the four metals studied, aluminium has the lowest 
specific heats in the range 0 to 20° K., and the electronic contribution forms 
an important part of the whole. Making use of the characteristic frequencies 
deduced in Part III from the specific heat data themselves as reported by 
Giaque and Meads and taking into account the electronic contribution, 
the specific heat curve has been computed. The graph thus obtained coin- 
cides with the observations reported by Kok and Keesom only at the lowest 
temperatures where the electronic heat contribution predominates. In the 
range from 5 to 20° K., the data reported by Kok and Keesom deviate pro- 
gressively from the theoretical curve, the observed values being everywhere 
in excess by quantities ranging from 8 to 20% of the calculated ones. They 
are also in excess of the values reported by Giaque and Meads for 15° and 
20°K. to a similar extent. 


Throughout the range now under discussion the specific heat of alu- 
minium is determined solely by the contribution of the residual spectrum 
with the addition of the electronic specific heat. The values reported by 
Kok and Keesom throughout the range would be in fair agreement with 
the theoretical calculations if the upper limit of frequency in the residual 
spectrum were taken as 146 cm.—! which is nearly the value deduced in Part II 
by an approximate method instead of as 160cm.-} calculated in Part III 
from the specific heat data themselves. But the former assumption would 
result in the theoretically calculated values at higher temperatures differing 
from those reported by Giaque and Meads by more than the probable 
magnitude of their experimental errors. Considering these facts, it appears 
correct to assume that the values reported by Kok and Keesom for alu- 
minium in the upper part of the temperature range between 0 to 20° K. are 
in excess of the real values by quantities of the order of 10%. 


SUMMARY 


The theoretical calculation of the specific heats of the four metals alu- 
minium, copper, silver and lead which in Part III were shown to be in 
satisfactory accord with the experimental data reported by Giaque and 
collaborators in the range from 15 to 300° K. have in the present paper been 
extended down to the absolute zero. When the contribution arising from the 
thermal agitation from the electrons is included, the calculated specific heats 
are in satisfactory accord with the results reported by various investigators 
on the low temperature specific heats of these metals. The agreement is 
excellent in the cases of copper, silver and lead, but not so satisfactory for 
aluminium, the available data for which are meagre and appear to be of 
doubtful accuracy. 
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ABSTRACT 


Many formule have been proposed to express numerically the rapid 
increase of the rotatory power of quartz with decrease of wavelength. 
An examination of these formule shows that they are wholly inappropri- 
ate for the case of quartz. It has been shown that the entire range of 
data from the visible to the extreme ultraviolet is accurately represented 
by a new type of formula involving only two constants, viz., 


ki 
P= 
where k = 7-19 and Ay = 0:0926283. A_ theoretical interpretation 
of the new formula has been given on the basis of a simple coupled oscil- 
lator model. 


INTRODUCTION 


THe phenomenon of optical activity was discovered in the year 1811 by 
Arago in pieces of quartz cut perpendicular to the optic axis. Soon after- 
wards, Biot showed that the rotatory power in the visible region of the 
spectrum varies inversely as the square of the wavelength, i.e., p = k/A*. 
The rotatory dispersion of quartz has subsequently been the subject of ex- 
tensive experimental study. This has been reviewed by Sosman in his book 
on the properties of silica (1927) and by Lowry in his treatise on optical 
rotatory power (1935). More recently, measurements have been pursued 
into the remote ultraviolet upto 1500 A by Servant (1939). It is evident 
from these determinations that Biot’s Law fails completely as we move away 
from the visible spectrum. Indeed, from 6000 A to 1500A the rotatory 
power increases by nearly 40 times instead of 16 times as required by the 
inverse square formula. 


There have been several attempts in the past at expressing the rotatory 
dispersion in terms of assumed characteristic frequencies. In the present 


* Now at the Crystallographic Laboratory, Cavendish Laboratory, Cambridge, England. 
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paper, these formule have been examined and it is shown that they are 
completely unsatisfactory for the case of quartz. A new type of formula 
has been suggested, a single term of which fits the whole range of data from 
the visible to the extreme ultraviolet very accurately. On the basis of a 
simple coupled oscillator model, a theoretical interpretation has been given 
to this formula. It is proved that the new formula is not open to the same 


theoretical objection which disqualifies all the earlier formule, which are 
principally of the Drude type. 


EARLIER FORMULZ FOR THE ROTATORY DISPERSION OF QUARTZ 


The first attempt at a theoretical formula for the rotatory dispersion of 


quartz was made by Drude in his Lehrbuch der Optik (1900), where he pro- 
posed the form 


p= 2 (1) 


where p is the rotatory power in degrees per millimetre, Q, is a constant 
corresponding to the characteristic absorption wavelength A, in microns, 
the summation being performed over all the r absorption frequencies. 
Drude found that the then available data could be fitted by a two term 
formula, where 


= 12-2, Q, = — 5-046, 
A,? = 0-010627, A,? = 0, 
the negative term being unambiguously indicated to be necessary. 
Lowry (1912) found that Drude’s formula is not accurate enough to 
fit his first series of extended measurements and hence altered it to 
_ 11-6064 43685 13-42 
— 0-010627 A? — 78-22” 


the additional third term corresponding to the infrared absorption band. 
The later and more accurate readings of Lowry and Coode-Adams (1927) 
required a further modification of (2), and their final formula is 


9563928113 
P = A®—0-000074 9° 1905. (3) 


where, instead of the single absorption band at 1031 A used in the previous 
formule, two new characteristic wavelengths at 1129A and 312A have 
been used, the effect of the infrared being substituted by a constant. 
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However, all these formule fail to represent the measurements of 
Duclaux and Jeantet (1926) extending up to 1854A. These authors have 
proposed the formula 

7-284 
P= 0-015114 (4) 


to approximately fit their data. But, as the discrepancy with this simple 
formula is quite considerable even in the near ultraviolet, another was 
worked out by Bradshaw and Livens (1929), viz., 


845 -694 0-40235 


™ (A? 0-01274912) (A? — 0-01274912)? 


838 -4320 0- 1331233 
(A? — 0:0120800) (A? — 0-0120800)? 


43-05794 2119-117 
+ 80-00) + — 80-00)? (5) 


It will be noticed that the rotatory power in this formula is expressed as a 
difference of two very large terms of nearly equal magnitude but of opposite 
signs involving two absorption wavelengths respectively at 1129 A and 1099 A 
and that the influence of the infrared has been retained. The formula, how- 
ever, deviates notably from the newer measurements of Servant (Joc. cit.) 
in the far ultraviolet region. Servant (1940) has therefore fitted his results 
with the expression 


+ 0-515 
(13-956 x 17639 x + 2-076) (6) 


p = 


where A is the wavelength of light in vacuum expressed in A. He has 
assumed the existence of two characteristic wavelengths, one at 1064 A, and 
another at 362A. The formula fits the entire data from 7000 A to 1500A 
well. Later, the same author (1941), realising that a multiterm formula 
involving numerous unknown constants is very artificial and arbitrary, has 
rejected (6) in favour of the simpler formula 


165-6 
p = (n,? + 2) — 1-325’ (7) 


where n, is the ordinary refractive index of quartz. The characteristic 
absorption assumed here lies at 1150A. The formula is approximately 
valid from 8000 A to 2000A and fails farther out in the ultraviolet. 
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Finally, Radhakrishnan (1947) has proposed the formula 


_ 4 4-617 2-311 
P= — 0-014161 * A2— 0-011236 — 0-000974 


0-000815 
— 0-1905. (8) 


This one is largely a refinement on the work of Lowry and Coode-Adams. 
The absorption at 312A is still retained, but the single term at 1129 A has 
been split into two terms at 1060 A and 1190 A. In addition a square term 
at about 1400 A, almost at the limit of continuous absorption of light in 
quartz, has been introduced. The formula gives a tolerable fit. 


Thus we see that a variety of formule involving numerous hypotheti- 
cal absorption frequencies have been suggested to express the rotatory 
dispersion of quartz. Even a cursory examination of these formule shows 
that most of them are extremely complicated and arbitrary in character and 
all, except one, are unsatisfactory in accounting for the actual observed 
facts. Apart from this, there is an even more serious objection which may 
be raised against them. 


KUHN’S SUMMATION RULE 


A fundamental theorem concerning the optical rotatory power of 
molecular systems was derived by Kuhn (1929) in the first instance by 
considering the behaviour of a model consisting of two coupled resonators, 
The theorem states that the sum of the numerators determining the contri- 
butions of the different absorption frequencies to the rotatory power should 
vanish, i.e., 


if 
then 
ar = 0 


Proofs of this summation rule involving no special assumptions regarding 
the molecular model have since been given by several authors (vide reviews 
by Condon, 1937; Kauzmann, Walter and Eyring, 1940; Mathieu, 1946). 


It is clear that the criterion furnished by this summation rule should be 
taken into account in deriving any formula, theoretical or empirical, which 


| 
‘ 
| 
| 


— 


The Optical Rotatory Dispersion of Quartz 151 


claims to express rotatory dispersion. Without entering into the detailed 
numerical calculations, even a rough estimate of 


Za (=F Q,/Ar”) 


in formule (1)-(8) shows that the summation rule is far from fulfilled by 
any one of them. As there is little justification in working out complicated 
expressions if we do not restrict ourselves to the basic requirement of the 
theory we may conclude that all the earlier formule are wholly inappropriate 
for the case of quartz. 


A New FORMULA 


The present writer has found (1952, 1954 4) that the rotatory dispersion 
of quartz from the visible to the extreme ultraviolet region of the spectrum 
is accurately represented by a simple formula involving only two constants, 
viZ., 

kA? 
The formula approximates to Biot’s Law when AS». Also for a given 
value of Ag, the rate of variation of the rotatory power with wavelength in 


this formula is much greater than in the Drude type of formula. The 
constants which the experimental data are 


k = 7°19, Ay? = 0-00858 or Ay = 0-0926283 p. 


(10) 


Here p is in degrees per millimetre and A is in microns. The correct proce- 
dure is to express A as the wavelength of light in vacuum. Servant has done 
so for his measurements in the remote ultraviolet, but the other authors have 
given the wavelengths in air. However, in the process of trying to fit a 
rotatory dispersion formula, Servant has effected the necessary corrections 
for a few wavelengths, extending from 1525A to 7000 A. The rotatory 
powers for these wavelengths along with the values calculated with the 
above formula are set out below (Table I). 


The agreement is very good from 6000A to 1525A. The error in 
this range, which never exceeds about two parts in a thousand, is not of a 
progressive nature. Data are available for many more wavelengths. The 
present writer (1954 c) has verified the validity of the formula for over 80 
wavelengths extending from 6000 A to 1520A. The fit is found to be ex- 
tremely good over this range, the error being sometimes positive, some- 
times negative. There is a slight discrepancy as we proceed towards the 
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TABLE I 
(iN p (expt.) p (calc.) p (expt.)—p (calc.) 
+1525 776-0 776-3 —0:3 
-1550 724-0 724-1 —0-1 
-1650 564-5 563-1 +1-4 
-1750 453-5 453-1 
-1850 374-0 374-2 —0-2 
-1950 315-5 315-3 +0-2 
-2000 291:3 291-3 +0:0 
+2250 205-9 205-9 +0-0 
-2500 154-60 154-56 +0-04 
-3000 97-72 97-61 +0-11 
-4000 50-26 50-17 +0-09 
-5000 30-86 30-84 +0-02 
-6000 20-93 20-96 —0-03 
-7000 15-13 15-20 —0-07 


infrared from the red end of the spectrum. This is no doubt due to the 
contribution of the infrared bands which have not been included in the 
above formula. Lowry’s measurements upto 3-7 definitely indicate the 
necessity of introducing a small infrared term. 


It can easily be shown that the single characteristic wavelength at 
0-0926283 » used in the above formula lies in the region of the spectrum in 
which the ultraviolet bands of quartz appear. A one-term dispersion formula 
for the ordinary refractive index of quartz has been fitted up using this same 
wavelength. The formula proposed is of the Sellmeier-Drude type which 
is theoretically valid irrespective of whether a polarisation field exists or 
not (Krishnan and Roy, 1956). The formula is 


. 
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where 
Ay? = 0-00858. 
A (in p) Ny (expt.) ny, (calc.) 
18547 1-676 1-673 
20006 1-649 1-649 
21107 1-634 1-635 
+2428 1-605 1-606 
1-577 1-578 
- 3587 1-564 1-564 
-4800 1-550 1-550 
-5893 1-544 1-544 
-7065 1-540 1-540 
*8447 1-538 1-538 


The agreement is reasonably good over the whole range. This indicates 
that the characteristic wavelength at 0-0926283, is approximately the 
effective average of the absorption spectrum of quartz in the ultraviolet. 
Of course, it is well known that the average characteristic wavelength opera- 
tive in refractive dispersion and in rotatory dispersion formule need not 
necessarily be the same, simply because the form of the expressions in the two 
cases are different (Condon, Joc. cit.), but in the particular case of quartz 
they happen to roughly coincide. 


The rotatory power normal to the optic axis was measured. for the first 
time by Voigt (1905). There have been a number of subsequent measure- 
ments, the most recent being those of Bruhat and Weil (1936), who have 
summarised their results as follows: the ratio of the rotatory powers of 
quartz perpendicular and parallel to the optic axis is independent of the 
wavelength for the spectral interval 5780 A- 2 40 A, the ratio being—0-54 
to within an accuracy of 0-5%. 


Evidently, the rotatory dispersion normal to the optic axis should also 
be expressible by a formula of the same type as (10) with A, unaltered and 
k reduced in the ratio of—0-54. 
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DIsCussION 


In their article in the Handbuch der Physik, Born and Goppert-Mayer 
(1933) arrive at the general formula 


2 Po 


It will be noticed that in addition to terms of the Drude type certain quadra- 
tic terms of the same form as suggested earlier are also present. We have 
now discovered that these terms play a predominant role in determining the 
rotatory dispersion of quartz. The simplicity of the new formula and the 
accuracy of its fit clearly shows that terms of the Drude type do not contri- 
bute in any observable measure to the rotatory power of quartz in the 
visible and ultraviolet region of the spectrum. 


After this work had been completed, it was found on going through 
some of the very early literature on the subject that in 1883 Lommel had 
suggested a formula of the type 


ar 
P= — 
for quartz and a few optically active liquids. He found that the available 
data for quartz from 8000 A to 2200 A were expressed tolerably well when 


log a@=0°8555912 or a=7-17119 and log A? = 7:9341257—10 or 
Ao = 0-0928032 


Lommel’s work has apparently been overlooked by most of the authors 
who have since been interested in this crystal. We have now found that this 
formula fits the data even up to the extreme ultraviolet. In addition, as we 
shall see later, we have been able to give it a theoretical basis while actually 


Lommel’s work was done before any proper theory of optical activity had 
yet been put forward. 


THE THEORY OF COUPLED OSCILLATORS 


The idea that optical activity arises from a particular spatial distribu- 
tion of coupled oscillators was developed by Born, Oseen, Gray, de Malle- 
mann and others (vide reviews by Condon, Joc. cit.; Kauzmann, et al., 
loc. cit.; Mathieu, loc. cit.). Later, Kuhn contributed greatly to the 
subject by putting forward a simple model which demonstrated all the 
essential features of the phenomenon. Kuhn’s model consists of two 
linear anisotropic resonators having natural frequencies of free vibration 
v, and vz respectively which are coupled with each other. Each resonator 
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is a charged particle which is assumed to be bound elastically to its own 
equilibrium position and capable of vibrating along a line. The analysis 
shows that the system has two modes of vibration, in each of which both 
resonators take part, one with frequency v, and the other with frequency v9. 
The contributions of the two modes to the rotatory power turn out to be of 
opposite signs when the frequency of the incident light is less than that of 
either resonator. 


In this model », and v, are assumed to be very different, so that the 
change of frequency brought about by the mutual coupling has been neglected 
altogether. This is, however, not permissible when v, and v, are very close 
to each other. In fact, in the limiting case when v, = v2, each frequency is 
actually split into a doublet, the magnitude of the splitting being dependent 
on the strength of the interaction between the resonators. 


We shall now discuss the optical behaviour of two identical linear 
anisotropic resonators coupled with each other. Let the natural frequency 
of vibration of the charged particle in each resonator be vp, when uncoupled. 
As a result of the coupling, v» would be split into two frequencies v, and v,. 
In one normal mode of vibration, the particles in the two resonators would 
vibrate with equal amplitude and in the same phase, and in the other normal 
mode they would vibrate with equal amplitude but opposite in phase. 

Let the first resonator be situated at the origin of the co-ordinate system 
and let its vibration direction be OX. Let the second resonator be at (0, 0, d) 
with its vibration direction parallel to OY. Suppose the charge and mass 
of the particle in each resonator to be g and M respectively. For convenience, 
we shall put q = fe and M = fm, where e and m are the charge and mass 
of the electron and f the oscillator strength. We shall call this coupled 
system the compound resonator. Let us examine the action of a light wave 
on such a compound resonator. 

Let right circular light given by 

x=E, cos wt 

sin owt 
be incident along OZ. The force exerted by the light wave on the charged 
particle of the first resonator is 


F, = fe Ey cos wt. 


The force on the charged particle of the second resonator which is rotated 
through 90° about the Z axis with respect to the first one is 


F, = — fe E, sin (wt — 94), 
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where 
2nitd 
fA being the mean refractive index of the medium. Since ¢ is small 
F, = — fe Ey (sin wt — ¢Cos wf). 


Now, the compound resonator has two normal modes of vibration and the 
force may be expressed in terms of the two normal co-ordinates as 


Re Fi + Fa 


Ry, = v2 (F, — F,) 
Substituting for F, and F,, we get 


+ 4) cos wt — sin wf] 


= a",, fe E, cos (wi +0), 


where (a",,)? = (1 + ¢), neglecting higher powers of ¢. Hence, the equa- 
tion of motion of the particles would be of the form 


+ cos (wt + 0). 


Therefore, 


fe Ey cos (wt + 0) 
(oF = 


Hence, the dipole moment induced 
= (a's, fe)* 


putting M= (fm. The dipole moment could be expressed as two compo- 
nents at right angles to each other, say mg and m,, where - 
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and 


These results are for a single compound resonator. If there are N such 
compound resonators all lying with their axes parallel to OZ but assuming 
different orientations about it, the mean polarisability per unit volume would 
be N (mg + mx)/2 Ep. 

Hence 

f 

(ny I)a, 2nN m w,? — 


In a similar manner, we get for the other normal mode 
where 
(a",,)? =1— 
Let us write 
np — 1 = (n,? — Ig, + (a? — 
For left circular light given by 
x = E,Cos wt 
y = Ep sin ot, 
we have similarly 


(nj? — 1)g, = (as, 


mM — w?’ 


(m? — 1)g, = (a’,,)? 


mM — 


where 
(ay)? = 1 — 4, 
(ay)? = 1+ 4, 
and 
my? — 1 = (my? — 1)g, + (my? — 
Therefore, 
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It will be noticed that the numerators in the two terms are equal in magni- 
tude but opposite in sign and this is in conformity with Kuhn’s summation 
rule. 


Let 
= wo? + 27% 
= wo? — 


We shall assume that the resonators are not very strongly coupled so that 
e is a small quantity. Putting 


and 


we get 


_ fe 


neglecting «? in the denominator. Hence 


p= — m) = 


where 


fe 


If the compound resonators were randomly oriented we would have to 
introduce a factor 4 in the expression for the rotatory power. 


If, as a more general case, we assume that the first resonator (0, 0, 0) 
has direction cosines a, 8, y and that the second resonator (0, 0, d) is rotated 
through an angle @ about the Z axis with respect to the first one, it may be 
shown that 

(13) 


where 


mNed sin 6 (a? + B?) A,tfe 
3 mc* 


It will be seen from (13) that the two enantiomorphous forms will have 
opposite rotatory powers, since pe< sin 9. 
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Thus, by means of an extremely elementary model we have been able 
to derive the correct expression for the rotatory dispersion of quartz. It 
should be exphasised that the new single term formula intrinsically satis- 
fies Kuhn’s summation rule and is therefore more justified than any of the 
earlier formule worked out for this crystal. 


We shall now try and interpret these results in terms of the structure 
of quartz. a-Quartz belongs to the space group D,* or the enantiomorphous 
D,°. The unit cell has three SiO, units spaced at equal intervals along the 
vertical screw axis, each of which is turned through 120° with respect to its 
predecessor. It is this screw structure which is responsible for its optical 
activity, the individual constituent units themselves being optically in- 
active. In view of the foregoing theoretical considerations, it can scarcely 
be doubted that the rotatory power of quartz arises primarily as a result of 
an interaction in the nature of a resonance between the similar polarisable 
units constituting the crystal. This type of interaction causes a splitting of 
the characteristic frequency of each individual unit and consequently, the 
rotatory dispersion formula is of the form (13). By a simple extension of 
the model which we have just considered, it may be shown in an approximate 
way (Chandrasekhar, 1954 c) that for an hypothetical crystal with a screw 
structure built up of such compound resonators, the rotatory dispersion 
formula for propagation of light along and normal to the optic axis are both 
of the quadratic form (13), a result which agrees with the experiments of 
Bruhat and Weil (loc. cit.). 


The same type of formula is also found to fit the data for cinnabar, 
benzil and sodium chlorate (Chandrasekhar, 1953, 1954 a, 1954 c) quite well 
and, in fact, much better than any of the previous formule proposed for 
these crystals. It may be remarked that, like quartz, the rotatory power of 
these crystals is present only in the crystalline state and not in solution. 


In conclusion, I have great pleasure in expressing my thanks to 
Professor Sir C. V. Raman, F.R.S., for the valuable discussions I had with 
him during the course of these investigations. I am also grateful to Dr. S. 
Ramaseshan for his kind advice and interest. 
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0. INTRODUCTION 


Tue classical study of the hydrodynamical behaviour of viscous liquids 
is based on the Navier Stokes Equations of motion, 


— grad p+pd4(u, v, w) (0.1) 
and the continuity equation 
div(u, v, w) =0. (0.2) 


The equations (0.1) follow by the application of Newton’s second law of 
motion to an element of the liquid and the use of the relation 


tj! = + Fod;' + Fid;! (0.3) 
between the stress tensor 7; and the tensor of rate of deformation 
bub 


where 8; are Kronecker deltas, 


ui, u®, u®; x}, x, x® stand for u, v, w; x, y, Zz 


respectively, and Fy, F, are functions of some material constants, 
as well as the invariants of the tensor d;'. Experimental evidence of the 
behaviour of some very viscous liquids noticed during the last decade led 
Reiner and others to examine the hydrodynamics of such general or non- 
Newtonian vicsous liquids. Reiner! has established thoroughly, that the 
most general relation between the stress tensor t; and the tensor of the rate 
of deformation d;‘ has the form 


t;! F,9;! F,dji 4+F od gid;*, (0.5) 
16] 
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which is a generalization of (0.4), including the additional term Fydqid;*. 
The presence of the coefficient F, gives rise to the phenomenon of ‘ cross. 
viscosity *. Writing Fy = — p, F, = 2p, F, = 2pe we recognize p, pw and 
He to be the mean pressure, the coefficient of viscosity and the coefficient 
of cross viscosity. 


1. THE EQUATIONS OF MOTION 


Adopting the relation (0.5) instead of (0.3) and replacing Fo, F,, F, 
in the above specified way, we write the equations of motion of non. 
Newtonian viscous liquids in the form 


by 
= p(X, Y, Z) — grad p+ (u, v, w) + 2pe div (dy'd;*). 


The vector div (d,'d;*) has the components 


t ve + 52) v, w) 

+ dads? + + + de2dy? + 


+ + (ds*)? + 
Substituting for d;' from (0.4) and using the continuity equation, viz., 
+ d,* + = 0, 


we obtain after simplification the following equations of motion in the 
Cartesian form: 
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dv 8 (du 
(du bw 
+ t+ + + 3) 
Su 


When the external forces belong to a conservative field ‘and are derivable 
from a force potential 2, we may write the equations also in the form 


+ wm 
+ (5 + 52) [oy + 5x 
+ 
we + ul 


+P + + vdv + [4 


bu 


dw , dy 5u\ 
bw dv\ 


i 
? 


dle 
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éw 
un + vé 


+ + 2) + Sy 


2. PLANE MOTION 


The above equations get remarkably simplified in the case of plane 
motion. If the motion is in planes parallel to the xoy plane we have w = 0 
and all the quantities are independent of z. The equations. of motion then 


assume the form 


— + + + dre (54) 


these equations can also be put in the form 


—(u, v, 0) 0, 


by 
. 
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From this we see by cross-differentiation, that the equation for vorticity in 
plane motions is given by 


(2.4) 


which is identical with the corresponding equation in the plane motions 


of classical or ordinary viscous liquids. It is significant that the cross-vis- 
cosity coefficient does not appear in this equation. 


if C= fq ‘ dr denotes the circulation of the fluid around any closed 
circuit, we may calculate d/dt C in the usual way. While this expression 
has an awkwardly long form in three dimensions, it assumes a simple, neat 
form in the case of plane flows. We have in fact, 


(circulation) = v J (dudx + Avdy) 


so that the circulation is conserved in plane motions if, and only if, 4udx 


+ Avdy is a perfect differential, or equivalently if the vorticity is a harmonic 
function. 


We may also examine here the conditions under which two non- 
Newtonian viscous plane motions are additive.* If and 
(Ue, Ve, fo, Po, $22) are the velocity components, vorticities, mean pressures 
and external force potentials respectively for two plane motions which are 
additive and = is the adjusted pressure, we have from (2.1), (2.2) 


Mylo 
— + 
— + + + 4%, { [8 


* After completion of the work in this paper we have noticed that A. C. Srivastava has 
also examined the question of ‘‘Superposability in Non-Newtonian Fluids,” vide Abstracts 
of the Second Congress on Theoretical and Applied Mechanics, New Delhi (October 1956). The 
details of his work are, however, not known to us. 
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Hence we obtain the following equation for determining the adjusted 
pressure 7. 


a(*) = + dx — + tals) dy — d (uy + V4V2) 


In order that this differential equation may lead to the determination of a 
single-valued function 7, + dx — (ul + dy must be an 
exact differential. Hence the condition for additivity of two plane motions 
is given by 


(Mabe + tabs) + § + vals) = 0, 


which is the same as in the case of the ordinary viscous liquids. 


3. We note here that the expressions arising in the equations of motion 
due to the cross-viscosity cannot be simplified any further than what we have 
recorded in (1.4), (1.5), (1.6) and as such the discussion of the general 
spatial flow with reference to vorticity, additivity of motion and circulation 
involves cumbersome, though straightforward analysis. However, on 
actually writing the equations determining the vorticity components we 
find that, unlike the plane motion, the cross-viscosity does make non-zero 


contributions to them. We are not reporting these equations on account 
of their bulk. 


4. DISSIPATION OF ENERGY 


If 5x Sy 5z denotes a rectangular parallelopipedal element of fluid at 
time ¢, with centre at (x, y, z) and t;' is the stress tensor, then the rate of work 
on the faces of the element by the tractions is given by 


u+ vy + t,!w) + 5, + t,2v + t,2w) 


+ (t,;3u + + | dx dy 


+ + + w dx dy 82 
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2 ow is 
ca 
= w, 5x dy + dx Sy 8z. 

w, 5x dy 5z gives the rate at which work is being done on the element in 
increasing the kinetic energy of the liquid while w,. 5x dy 5z measures the Vi 
rate at which work is being done in the deformation of the element. It 
Writing fe 
= — pd; + + 
( 
we see that 
We = (— p + 2udi! + + (— p + 
27a) OV 3 OW ‘ 
+ 4°) 5 + (—P + 2uds + d3*) 
( 
+ (2pds? + 2ped (2d,") + (2Qud,* + 2peda*d;*) (2d,*) 
+ (2uds! + (2ds"), 
where the dummy index a, as usual indicates summation. Since F 


+ dt + 4 = 0, 


we may simplify w. and finally put it in the form 
We = {2 (dy')? + 2 (dy)? + 2 (ds)? + + (2d,%)? + 
+ Ope + — (ds*)* — (ds)? 
— (d,})?}. 


The coefficient of 64¢ can be written in the following determinental form. 


In the case of plane motion we see again that the cross-viscosity disappears 
in the expression for w, for d,! = d;* = = 0. The above determinant 


1 1 1 
ds ds 
2 2 2 3 
- 
a, ds 
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is evidently not of a definite sign, unlike the coefficient of » which is a positive 
definite quadratic form in the elements d;'. Since by its physical nature ws 
cannot be negative, we must have 


6 | He + — dy" (ds*)* — (ds!)* — 
<p {2 (dy)? + 2 + 2 + (2d5*)® + + 


It is worthy of note here that the coefficients of viscosity (4) and cross- 
viscosity (ue) are not on equal footing, in the analytical or physical sense. 
In the first place they do not have the same physical dimensions. The 
features of viscosity and cross-viscosity are indeed quite distinct. While 
the coefficient of » in the dissipation term is a positive definite quadratic 
form in the tensor components d;', the coefficient of 4e in the same term is 
a cubic form in the elements d;' with no definite sign. 


The term F,d,‘d;* involving the cross-viscosity 4e = $F, has been intro- 
duced in the stress-rate of deformation relation (0.5) to explain the hydro- 
dynamical behaviour of highly viscous liquids and the cross-viscosity term 
(F.d,'d;*) in (0.5) may be regarded as a correction to the classical relation 
(0.3) in the case of ordinary viscous liquids. The dissipation term w, should 
of course always be positive, but this need not be the case for the part of 
W, arising due to cross-viscosity. This forces the conclusion that we cannot 
conceive of a viscous fluid motion in which the cross-viscosity terms alone 
are important and the ordinary viscosity terms are neglected. 


5. Several interesting questions have been examined in the flow of viscous 
incompressible fluids by the so-called ‘ inverse ’ and ‘ semi-inverse ’ methods. 
An investigation of a partial differential equation is called ‘inverse’ if 
the boundary conditions in it are not prescribed at the outset and additiona] 
properties of the field (of analytical or geometrical or physical nature) are 
invoked. An investigation is semi-inverse if some of the boundary condi- 
tions are prescribed in the beginning, while others are left open, only to be 
obtained indirectly from simplifying assumptions about the properties of the 
field. 


Some of the major investigations in the case of viscous flows employing 
the inverse or semi-inverse method deal with problems of a characterizational 
nature. Thus Jeffery? proved that the only steady plane rotational flows 
of a viscous liquid (in the absence of body forces) for which the isocurls 
(lines of constant vorticity) coincide with the stream lines of a plane irrota- 
tional inviscid flow are those in which only concentric circles or parallel 
straight lines are admissible as stream lines. The investigations of Jeffery? 
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and Hamel*® have also shown that the only steady, plane rotational flows 
of viscous liquids (in the absence of body forces) for which the stream lines 
coincide with the stream lines of a plane, irrotational inviscid flow are the 
spiral flows (or their degenerate cases). Problems of this nature require 


the use of the equations wherein the terms involving pressure have been 
eliminated. 


If such investigations (are to be extended) to the case of a non-Newto- 


nian viscous liquid, one has to start with the equations (2.3). If w(x, y, 1) 
denotes the stream 


function so that the velocity components are 


“w= 


(2.3) can be put in the form 


5x2 dy? 


(5.1) 


On eliminating the pressure term by cross-differentiation we notice that the 


cross-viscosity term also drops out as is evident from (2.4) and we have 
merely 


5 Ap) _ 
57 + vA Ad. (5.2) 


The first term will not be present in a steady flow. Since this equation does 
not contain the cross-viscosity coefficient and this is the starting point for all 
inverse and semi-inverse investigations, we may observe in passing, that such 
investigations with viscous steady plane flows will not suffer any modifica- 
tions in the case of non-Newtonian fluids. 


5. SUMMARY 


The paper deals with some general considerations in the motion of 
highly viscous liquids for which the stress-rate of deformation relation is 
t;' = F,8;' + F,d;' + F.d,'d;*. The equations of motion, the equations 
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for vorticity and the energy dissipation term are obtained here for a general 
spatial flow. In the special case of the plane flow, the vorticity equation is 
free from the cross-viscosity coefficient. This simplifies the study of 
problems like the addition of two plane motions, consideration of all 
circulation preserving motions as well as some ‘ semi-inverse ’ investigations 
in plane motions of non-Newtonian liquids. 
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QuiTE a number of interesting relations have been noticed in recent years 
(cf. the reference at the end of this paper) concerning the ultraspherical 
(and other orthogonal) polynomials, in connection with the proof of 
Turan’s inequality which had an extraordinary resonance in the recent 
literature on Special Functions. In the course of proving Turan’s inequality 
for Legendre polynomials, viz., 


An (x) = [Pn — Pry (x) Pri 0, @> 1, |x|< (1) 


B. S. Madhava Rao and V. R. Thiruvenkatachar! established the convexity 
property of 4y(x) by showing that 


2 

V. R. Thiruvenkatachar and T. S. Nanjundiah? employed the formula 

(1 — x’) D,»™ (x) = kn,y An,» (x) (3) 
where 

P, (x) 2 (x) (x) 

An, A(x) = ( 5) (5 (Gant (1) (4) 

(x) = Wo) = (x) ex); (5) 
and 


Kn, = n(n + 2A) 


for proving Turan’s inequality for ultraspherical polynomials P,,™ (x) in 
the form 


4n x(x) 20, AA>—4, nol, |x|<1). (6) 


In his work on explicit evaluations of Turan expressions, A. E. Danese? 
has obtained a number of relations which substantially contain the proof 
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of Turan’s inequality for the concerned orthogonal polynomials. K. Venkata- 
chaliengar and the present author have shown that if 


A,™ (x) = [P — P na? (x) (x), 


then “4 P, (x) can be represented as a numerical multiple of the 
Wronskian of Pyyy (x) and Py, (x) and have also obtained a series 
representation of positive functions (apart from a multiplier) for = An™ (x). 
The aim of the present note is to point out that these several relations 


employed or encountered in the course of proving Turan’s inequality for 
P, (x), are characteristic for the ultraspherical polynomials. 


THEOREM 1|.—If 
An (x) = [Pn — Pris (x) (7) 
then 


ad 2 d 2 
dx? 4n (x) = — n(n +1) | (8) 


This relation is characteristic for Legendre polynomials, i.e., If{ fy (x)} 
is a sequence of polynomials such that the degree of 


and 
{L fn CO)? — Sra O)} = Ex (9) 
then 


fn (x) = Pn(x), (n=0, 1, 2,....). 
Proof.—We notice that 

So (x) = Po (x), (x) = Pi (2). 
Let us assume that 

Sr (x) = Py (x), 
and let 


Comparing the coefficients of x?"-* on both sides of (9), we find that 


= I(n >) (n +- 1)! 
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the highest coefficient in Py, (x), so that we may write 
Fria = (x) + (x) 
where g(x) is a polynomial of degree <n. This substitution for Sra (x) 
in (9) leads to the equation 
dx? [g (x) Pana (x)] = 9, (n = 1, 2, 3,....). 


This second order differential equation has the solution 
& (X) = + (x), (n 2 1) 


which is clearly not compatible with the polynomial nature of g(x). We 
conclude that c; = c,=0 or g(x) =0 and thus prove the theorem. 


THEOREM 2.—If 


4, (x) = [P,™ — (x) (x) (10) 
then 
d Pria™ (x) Pra” (x) 
4,™ (x) = 
dx 


(x) dx (x) 


and the relation is characteristic for the ultraspherical polynomials P, (x), 
If {fn (x)} is a sequence of polynomials such that the degree of 


= 03 (X) = 1, AM = 2Ax 


and 

then 


fa™ = Px™ @), @ = 9, 1, 2, ....). 


Proof.—We take 4 > — 4 but omit the case A=1 for then we _ have 
4, (x) = 1 and both sides of (11) are zero. As in the proof of Theorem 1, 
if we take 


(x) = (x), O<r<n 


fro 


é 


and 
| whe 

We 
tha 
7 ; 


11 (x) 


12) 


ve 
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and set 
by comparing the highest coefficients in (12), we notice that 


(*) = Pra (x) + (2), 


where g (x) is a polynomial of degree <n. Employing this relation in 
(12) and using (11), we arrive at the relation 


(n+ 1)g™ (x 4 (x) 


(A) 
Op, =0, (13) 
We may evaluate g(x) by direct integration. Otherwise, remembering 
that Py» (x) and d/dx Py_,™ (x) have no common factor, we conclude 
from (13) that g™ (x) = aPy_,.™ (x). It follows then that 
a[@ + 1) + @+ 2A— Pra™ (x) = 0. 


Thus we are forced to accept that a = 0 and the theorem is proved. 
THEOREM 3.—The difference equation 


+ 2A — 
(x) = (n+ 2 (x) 


and the series representation 


d 
+ 1)/2] 
(n — 2« + 2e +A4+1) 
In — +2A+1 
X [Pn (15) 


are characteristic relations for the ultraspherical polynomials P,™ (x). The 
proof runs exactly on the same lines as in Theorem 2. Here also, when 
A=1 both the members of the equation reduce to zero, 


We 
(x), 
_| 
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THEOREM 4.—The formula 


(1 — x’) iF: (x) 4 Pr 


= n(n +- 2d) (x)]? — (n + 1) (n + 2A — 1) 
x (x) (x) (16) 


is characteristic for the ultraspherical polynomials. Here again the proof 
is as in that of Theorem 2. In the same way we may show that the relation 


(mn + 1) + 2A — 1) — (x) (x)} 


+ (2A — 1) [Pr™ (x)]’, 
where 


is characteristic for the ultraspherical polynomials. 
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ABSTRACT 


In the present paper we have established uniqueness of maximum 
pressure for composite charges, when each component charge has the 
cubic form-function. The conditions for the occurrence of maximum 
pressure in or at the end of any stage have been explicitly obtained. 


1. INTRODUCTION 


To determine the stage or stages of burning in which pressure-maxima can 
occur and in particular to determine whether the maximum pressure would 
be unique are important problems from the point of view of designing of 
guns for use with composite charges. Kapur (19564, 5) established the 
uniqueness of maximum pressure under the following assumptions: 


(i) rate of burning proportional to pressure for each component charge. 
(ii) = Ye = = 
(iii) Specific volume = Co-volume for each component charge. 
(iv) the standard quadratic form-function of the type 
z=(1—f)(1+ f) (1) 
for each component charge. 


The first two are the usual assumptions. The effect of taking co- 
volume-correction terms has been examined in a recent paper [Kapur (1957)], 
where it has been proved that except for very highly progressive forms 
obtained by inhibiting certain surfaces from burning—forms which are 
used in rockets but not in guns—the uniqueness of maximum pressure 
holds even when we take these co-volume-correction terms into account. 
As for assumption (iv), most of the forms used in practice can be more 
correctly expressed by the cubic form-function 


z=(1—-f)( + (2) 
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Since in most of the common forms, tube, cord, ribbon, flake, y is 
small and positive, form (1) is preferred for analytical solutions, specially 
because form (2) gives rise to complicated solutions. It is however desir. 
able to establish the uniqueness theorem under as general conditions as 
possible and fortunately we find, that the proof, while deeper, is not more 
complicated for the cubic form-function. We have also to remember that 
for highly degressive shapes like spheres % can be as high as unity. 4% can 
also be negative for inhibited charges, but such shapes are not used in guns 
and consequently we shall assume % to be positive in our discussions. The 
method used here to establish the uniqueness of maximum pressure for 
the cubic form-function is the Equivalent Charge Method of Corner (1950) 
and Clemmow (1951) as improved by Kapur (1956 c). 


2. THE Basic EQUATIONS FOR THE EQUIVALENT CHARGE 


Let F;, Cj, Dj, Bi, 24,f;, 9, 4; refer to the ith component charge, and 
let F, C, D, B, z, f, 6, & refer to the equivalent charge, then the fundamental 
equations for the equivalent charge [Kapur (1956 c, 1957)] are: 


FCz = Ap(x + D+ 4(y — 1) wy? (3) 

Da, = — Bp. (4) 
42 = w, 9? = Ap (5 
z= A, + —f) — —f)? + Dy (1 — Sf) (6) 

when 

1 1 

where 
Ar= (7a) 
By = Aki (1 + + Ha, (7) 
E, = (0; + 21), (70) 


i=r 


(74) 


i=r 


% 
¢ 
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is and 

ally FC; 

as z 

ore 

hat _ _ Bil Di Bil Di B) 

pr = "BID = 

uns 

The and charges have been so numbered that 

for 

50) B,’ > B,' > (10) 
Using (4) and (6), we get 

nd 

v= Bw, (11) 


where f) is value of f at shot-start. 


From (3), (4), (5), (6) and (11), we get 


4) A(x+D dp 
dp Br + —f) — — 


+yM(fo—f), (12) 
6) where 
A2D?2 
M = BFC (13) 
a) 


is the central ballistic parameter for the equivalent charge. 


Also from the equations of the rate of burning 


| D, df__ 
dt B dt ? 

bod 

c) 

or 
d) (1—fy) (14) 
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Using (14), we can write (12) as 


A(x +1) dp +f 6D, 
“FO df= ft + (yM + 2, 
2E 3D, yM 
(15) 
3. Discontinuity OF dp/df AT THE END OF A STAGE 
At the end of the r-th stage, f= 1 — 1/k, and therefore from (12) 


dp _ FC _2E, 3D, yM 


(16) 
where x, is the shot-travel till the end of the r-th stage. 


Similarly at the beginning of the (r + 1)th stage 
Brus » 


df k k,* 
+7 — — (17) 
therefore increment in dp/df at the end of the r-th stage 


r 2 3 
= A(x (Bria — Br) + (Eri. — — 


(Dry — D,)| 


Now if 
&ft+ (19) 


is to represent a form-function for a charge, dz/df should be negative for 
all values of f lying between 0 and 1, though it may vanish at f= 0 and 
f= 1. From this fact, we easily deduce 


(ii) 1 + 6+ 4,>0. (20 b) 


From (18) and (20a), we find that in crossing from the r-th stage to 
the (r + 1)th, dp/df always receives a non-negative increment, Therefore 


if 
q 
t 
(i) 1— 6,>0 (20a) 


1) 


5 
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if dp/df is negative at the end of a stage, it can become positive at the begin- 
ning of the next; but if dp/df is positive at the end of the previous stage, it 
cannot become negative at the beginning of the next. 
4. UNIQUENESS OF MAXIMUM PRESSURE 
We shall establish the uniqueness for the case when 
tn >0 (21 a) 
, On> vn (21 


These conditions are satisfied by most charges used in practice, as is 
shown by the following table: 


Charge Cord Tube Sphere Cube Ribbon Square-flake 
1/A 0 1 1 


Here A and ware large and are above 100 for service propellants. 


For this case 6, > 0 for all r, and from (7), A,, B,, E,;, D,, are positive 
for all r. Also 


= yM +2 3 dike | + — 


i=r 


Since 
ki< ky for r>n, 


we find from (21) and (22) that 


ae > 0 for all r. 
ky 
Now let 
K, =B, — 2Er + yM (1 — fo), (23) 
ke? kr 


then (15) shows that 


(i) If K, > 0, dp/df does not vanish in the r-th stage, and pressure 
continues to increase steadily in this stage. 


6) 
6D, 
yM + — 
ky 
8) 
9) 
or 
nd 
a) 
6) 
re 
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(ii) If K, = 0, dp/df can vanish for a negative value of f, (which is in- 
admissible) and for /, = 0 (i.e., at the end of the r-th stage). 


(iii) If K, < 0, dp/df can vanish for two values of f,. One of these 


is negative and therefore inadmissible. The other is positive and would be 
admissible if it is less than unity. 


Thus dp/df can vanish at most once during a stage and a pressure 
maximum or minimum can arise at this instant. However (15) shows that 
since D; and yM + 2E, — 6D,/k, are positive «nd since f, has some posi- 
tive value at the beginning of the r-th stage and is zero at its end, the 
negative contribution to dp/df is reduced during the r-th stage, and there- 
fore if dp/df is negative at the beginning of the r-th stage, it can become 


positive at its end and not vice-versa. Consequently only a pressure maxi- 
mum can arise during a stage. 


At the end of a stage also, since dp/df receives a non-negative increment, 
it can change there from negative to positive and not vice-versa, and thus 


a pressure maximum can occur there, but a pressure minimum cannot occur 
there. 


Thus we find that both during a stage and at the end of a stage, dp/df 
can change only from negative to positive (or zero) and consequently once 
dp|df has become positive, it cannot become negative again. The unique- 
ness of maximum pressure is thus established. 


5. CONDITIONS FOR THE OCCURRENCE OF MAXIMUM PRESSURE 
IN OR AT THE END OF ANY STAGE 


The maximum pressure would agcur in the r-th stage if dp/df is nega- 


tive at the beginning of this stage and positive at its end. The correspond- 
ing conditions are, from (12): 


Ky k 


The maximum pressure would occur at the end of the r-th stage if 


(i) dp/df is negative at the end of the r-th stage and positive 
beginning of the (r+ 1)th, or 


at the 


(ii) dp/df is zero at the end of the r-th stage, or 
(iii) dp/df is zero at the beginning of the (r + 1)th stage. 


i 
} 
; 
: 
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The corresponding conditions are: 


and 


3D M 
— Bra + — yM (1 — fo) > 0 


(ii) —B, + — fo) =9 
ky ky 


(iii) —Brat 3Dris yM (1 — fo) = 9 


The maximum pressure would occur at all-burnt if 


Remembering that - = 1 and using (7), this simplifies to 
An (1 — On) > yMfo. (26) 


When #, = f. = .... = tn = 0, so that D, = 0 for all r, the condi- 
tions (24), (25) and (26) reduce to the corresponding conditions for the 
standard quadratic form-function [Kapur (1956 d)]. 
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INTRODUCTION Di 
Jessop! has reported that in alkaline triethanolamine solutions nickel gives 
a wave of poor shape. Preliminary investigations in this iaboratory indi- 
cated that the use of ethanolamines under varying conditions produced Tr 


waves which could be of great help in the polarographic estimation of nickel. 
Consequently, a detailed investigation of the polarography of nickel was 
undertaken and the results obtained are reported below. 


EXPERIMENTAL 


The polarograms have been recorded automatically by the Cambridge 
Instrument Company polarograph. The procedure and other experimental 
details are substantially the same as those described previously.*, The drop- 
time of the capillary has been 5 seconds in | N potassium chloride and the 


value of m= 1-150 mg. per second. The temperature of recording the 
polarograms has been 30° + 0°1°C. 


RESULTS 
1. Effect of pH 


In these experiments, the concentration of the amine was kept at 0-1 M, 
and that of nickel at 1-6 millimolar. Gelatin (0-005%) was employed to 
suppress the maxima which appear in the solutions containing mono- and 
diethanolamines. Since the solubility of the complex in 0-1 M solutions is 
low, slight precipitation occurs at about pH 10. It is surprising that at pH 11, 
both in mono- and in diethanolamines, maxima persist even after the addition 
of gelatin under the above conditions. Table I presents the results of ana- 
lysis of the waves. 
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TABLE I 


Effect of pH on the polarographic behaviour of nickel in presence 
of ethanolamines 


Half-wave potential 
Amine used E, vs. S.C.E. 
(Volt) 


Monoethanolamine. . —1-038 
—0-862 
Maximum not suppressed 


—1-043 
Diethanolamine .. —1-069 
—0-933 
Maximum not suppressed 
—1-163 
Split waves 
Triethanolamine  .. —1-093 


A perusal of the table indicates that the reduction process is irreversible 
in most cases except in alkaline monoethanolamine solutions. The pheno- 
menon of “ split-wave’’ which is absent in monoethanolamine, becomes 
noticeable only in highly alkaline diethanolamine solutions whereas in tri- 
ethanolamine solutions it starts from a pH 7-20. The complexity of the 
wave in alkaline triethanolamine solution can be seen from Fig. 1. 


2. Effect of Variation of the Composition of the Base Electrolyte 


During the present studies noteworthy features of complex formation 
have been revealed only in alkaline region as in the case of cadmium.* The 
results obtained with solutions of sodium hydroxide, sodium carbonate and 
ammonium chloride-ammonium hydroxide are presented below. 


(i) Sodium hydroxide.—The polarograms in Figs. 1 and 2 show the effect 
of sodium hydroxide on the behaviour of nickel in presence of ethanolamines, 


185. 
—E; 
(Volt) 
0-108 
0-042 
12-70 —0-986 0-029 
0-029 
0-116 . 
| 0-048 
0-099 
7:20 —1-100 0-087 - 
8-65 —1:065 0-076 
—1-196 0-038 
9-70 —1-090 0-065 : 
—1-207 0-034 
12-70 ~1-240 0-074 
—1-386 
‘13-50 —1-292 0-068 
—1-490 
~1-680 
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where 0-005% gelatin has been employed to suppress the maxima in solu. 
tions containing mono- and dietiianolamines. 
ethanaolamine (Table ID) it is evident that a single reduction wave gives 


TABLE II 


Effect of sodium hydroxide on the polarographic behaviour of 
nickel in presence of ethanolamines 
Concentration of nickel: 1:60 millimolar 


As we go from mono- to tri- 


Amine used 


Concentration of 


Sodium 


Half-wave 

potential 

Amine hydro- vs. S.C.E. 
id. 


xide 
(Molar) (Molar) 


(Volt) 


E; —E; 


(Volt) 


Wave height 
at 5/30 


Monoethanolamine 0-50 


Diethanolamine .. 0-30 


Triethanolamine .. 0-10 


1-00 


1-00 


=? 
—1- 
wale 
ons 
wits 


0-030 
0-028 
0-033 
0-034 


0-048 
0-050 
0-045 
Irreversible 
waves. Mea- 
surement 
difficult 


0-074 
0-048 
0-057 
0-063 
0-068 


0-061 


0-071 


0-068 


45-1* 
46-3 
38-7 

40-8 
41-7 
40-9 
39-1 
35-3* 
36-6 
34-1 
37-5* 
40-7 
36:8 
34-2 


Measurement 
difficult 


* Values are approximate. 


sl 
re 
re 
i 
1 
28 
1-00 )45 
0-50 84 
1-00 102 
97 
0-65 125 
1-00 138 
0-30 150 
10 
0-65 7 270 
| 80 
1-00 1-00 B80 
90 
390 
0-30 0-10 280 
420 
0-60 0-10 320 
500 
1-00 0-10 330 
510 
0-10 1-00 292 
500 
: 670 
0-30 1-00 333 
540 
690 
0-60 1-00 365 
530 
740 
550 
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place to ill-separated reduction steps of two or three different species in 
sluggish equilibrium.* It is clear from the table that reduction proceeds 
reversibly in monoethanolamine, while in di- and triethanolamines, the 
reduction process is irreversible. These results are in contrast to the results 
obtained with cadmium,* where the complexes of diethanolamine and of 
monoethanolamine showed greater similarity. 


(ii) Sodium carbonate.—In Fig. 3 are given polarograms of nickel in 
monoethanolamine. Curves in di- and triethanolamines are similar to those 
in monoethanolamine. Gelatin (0-005%) is added to the base solutions 
containing mono- and diethanolamines to act as the maximum suppressor. 
Table III gives the half-wave potentials and the values of E; — E, for the 
systems studied. 


The results given in Table III indicate that the reduction of the mono- 
ethanolamine complex proceeds very nearly reversibly while in di- and tri- 


TABLE III 


Effect of sodium carbonate on the polarographic behaviour of 
nickel in presence of ethanolamines 
Concentration of nickel: 1-60 millimolar. 


Concentration of 
Amine used Half-wave 


Sodium potential E;—E;, Height of the 
Amine carbo- E, vs. S.C.E. wave s/50 


nate 
(Molar) (Molar) (Volt) (Volt) 


Monoethanolamine 0-034 
0-032 
0-033 
0-034 
0-032 
0-033 


0-068 
0-064 
0-068 
0-076 
0-069 
0-064 


ooo 


Diethanolamine 


eeesss sss 


SSRSER 


Triethanolamine .. 


0-095 Precipitation 
42- 


0-088 2°5 
0-068 34-2 


SsS 
COO HHH OOO HHH 


ooo 


solu- 
) tri- 
: 
: 
: 
‘ 
i 
I 
nt 
—1-190 
= —1-290 
A4 
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ethanolamine solutions the reduction process is irreversible in character. 
Even in sodium carbonate solutions, the nickel-diethanolamine complex 


resembles the triethanolamine complex, thus differing substantially from 
cadmium’ complexes. 


(iii) Ammonium chloride (1 M) and ammonium hydroxide (1 M).—The 


effect of ammonium chloride-ammonium hydroxide mixture is given in 
Table IV. 


TABLE IV 


Effect of ammonium chloride (1 M)—ammonium hydroxide (1 M) on 
the polarographic behaviour of nickel in presence of ethanolamines 


Concentration of nickel: I-60 millimolar. 


Concentra- Half-wave 
Amine used tion of potential E; —E; Height of 
amine E; vs. S.C.E. the wave at 
(Molar) (Volt) (Volt) s/50 


Monoethanolamine .. 


—1-066 0-059 63-2 
0-50 —1-081 0-056 62-4 
—1-226 
1-00 —1-092 0-052 59-8 


—1-275 


Diethanolamine 0-10 —1-070 62-3 
0-50 —1-097 0-063 58-9 
1-00 —1-114 0-064 55°5 
Triethanolamine ue 0-10 —1-080 0-062 57°7 
0-30 —1-132 0-077 55-9 
0-60 —1-173 0-098 52-4 
1-00 —1-189 0- 47-9 


It is clear that the reduction proceeds irreversibly in all the cases. With 
mono- and diethanolamines the irreversibility is almost constant between 
0-1 and 1-0M concentration of the amine, while in the case of triethanol- 
amine it increases with increasing concentration of the amine and becomes 


almost constant when the concentration of the amine is increased beyond 
0-6 M. 


In the case of monoethanolamine another small wave appears at — 1-23 V 
vs. S.C.E. in 0-5 M solution. This value increases to ca. — 1-28 V vs. S.C.E. 
in 1M solution. With diethanolamine, however, the bend before the diffu- 


(2) 


(3 


sion 
hand, 
in 
i 
of nic 
_ 
(1) . 
( 


Polarographic Behaviour of Metals in Ethanolamines—III 189 


sion current plateau is pressed inside. Triethanolamine complex, on the other 

hand, does not show such phenomenon. 

3. Effect of Concentration of Nickel on the Diffusion Current Constant 
in Presence of Sodium Hydroxide, Sodium Carbonate and Ammonium 
Chloride-Ammonium Hydroxide 
The results obtained are given in Table V. The effect of concentration 

of nickel on the diffusion current constant has not been tested in the case 

TABLE V 


Effect of concentration of nickel on the diffusion current constant* 
Concentration of the amine: 1 M 


Concentration Diffusion Concentration Diffusion 
of nickel current of nickel current 
(Millimolar) constant (Millimolar) constant 


la 
tt 


(1) B.S.: 1-00 M Moen+0-10 M NaOH (5) B.S.: 1:00 M Moen+0-10 M 
+ 0-025% gelatin Na,CO, 
X = 1-415 at — 1-30 V vs. S.C.E. + 0:025% gelatin 
1-600 2-83 X=1-424 at —1-30 V vs. S.C.E. 
8-000 2-80 1-600 2-99 
(2) B.S.: 1-00 M Dien+0-:10 M NaOH 8-000 3-01 
+ 0-025% gelatin 16-00 2-93 
X = 1-408 at — 1:40 V vs. S.C.E. (6) 1-00M Moen+1-00M 
1-600 2: 37 Na,CO, 
8-000 2-35 + 0-025% gelatin 
(3) B.S.: 1-00 M Dien+1-00 M NaOH X= 1-424 at — 1.30 V vs. S.C.E. 
+0-005% gelatin 1-600 2-34 
X = 1-408 at — 1-40 V vs. S.C.E. 8-000 2-32 
1-600 2-32 16-000 2-32 
8-000 2-30 (7) B.S.: 1-00 M Dien+-0-10 M Na,CO, 
(4) B.S.: 1-00 M Trien+0-10 M NaOH + 0-025% gelatin 
X = 1-373 at — 1-60 V vs. S.C.E. X = 1-424 at — 1-30 V vs. S.C.E. 
1-600 2:46 1-600 2-16 
3-200 2°44 8-000 2-14 
8-000 2°44 (8) B.S.: 1-00 M Dien+1-00 M Na,CO, 
+ 0-005% gelatin 
X = 1-424 at — 1-30 V vs. S.C.E. 
1-600 2-04 
8-000 2-03 
(9) B.S.: 1:00 M Trien+0-:10 M 
X = 1-383 at — 1-60 V vs. S.C.E. 
1-600 2°44 
3-200 2-46 
8-000 2-43 


*Abbreviations :-B.S.—Base solution ; X-ni in mg! sect Moen-monoethanolamine ; Dien- 
diethanolamine ; Trien-triethanolamine. 
AB 
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of ammonium chloride and ammorium hydroxide since nickel is quite soluble 
in ammonium hydroxide-ammonium chloride and the incorporation of 
ethanolamines does not in any way help in the practical estimation of the 
metal. In solutions containing sodium hydroxide, | molar amine is incorpo- 
rated in the base solution since the nickel complex is not sufficiently soluble 
at lower concentrations of the amine. Since in certain base solutions maxi- 
mum appeared at higher concentration of nickel when 0-005% gelatin is used, 
the concentration of gelatin has been increased to 0:025%, though this brings 
about a decrease of the diffusion current by nearly 12%. In triethanolamine 
solutions, normal sodium hydroxide gives rise to waves of poor shape, while 
in the case of monoethanolamine solutions the nickel gets precipitated due to 
the low solubility of the complex. Hence it can be concluded that the use of 
1 M sodium hydroxide in base solutions is not suitable for analytical estima- 
tion of the metal. As in the case of sodium hydroxide 0-025% gelatin had to 
be employed to prevent the appearance of maxima at higher concentration 
of nickel in some of the base solutions. From the results it can be concluded 
that in sodium carbonate solutions also, nickel can be estimated in all the 
three ethanolamines. The use of monoethanolamine has a decided advantage 
over di- and triethanolamines as the reduction process is reversible and the 
value of diffusion current constant is greater than in the other two amines. 


DISCUSSION 
1. Splitting of the Polarographic Wave in Triethanolamine 


The splitting up of the polarographic wave is quite prominent in sodium 
hydroxide solutions (Fig. 1). The main wave in each of the polarograms is 
due to the complex containing the amine and hydroxyl while the other two 
are due to complexes containing more of the hydroxyl in sluggish equili- 
brium, analogous to the complexes of cadmium® studied by the present 
author. Of the ill-separated doublets that appear between pH 10 and 7, the 
first corresponds to the reduction of aquonickel ion as its half-wave potential 


is found to be — 1:06 V and the second corresponds to a complex of the 
amine (Table J). 


2. Splitting of the Polarographic Wave in Mono- and Diethanolamines 


In contrast to the behaviour of triethanolamine, splitting is practically 
absent in monoethanolamine solutions. However, the small second wave 
that appears with ammonium chloride-ammonium hydroxide solutions is 
presumably due to a species of the complex, containing more of the amine 
than of the complex giving rise to the main wave. Similar explanations can 
also be given for the elongated bend with the diethanolamine solutions. 
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The diethanolamine solutions show the splitting in 1-0 M sodium hydroxide 
solutions. In this case also it is reasonable to expect the second species, to 
contain more of the hydroxyl groups. 


3. Effect of pH 


Mono-and diethanolamines.—In acid solutions (pH ca. 6) the half-wave 
potential is of the order of — 1-04 V vs. S.C.E. and the wave is highly irre- 
versible the value of E,;— E, being of the order of 0-100 V (vide Table I) 
as compared to a value of 0-029 V required by a reversible two-electron reduc- 
tion process. The half-wave potential of hexa-aquonickel ion is ca, — 1-10 V 
ys. S.C.E. and the wave is irreversible. It therefore appears that in acid solu- 
tions it is the aquonickel ion that is getting reduced. When the pH is 
increased, the half-wave potential and the value of E;—E,; drops. This 
indicates complex formation as it is known‘ that nickel is more reversibly 
reduced in presence of complexing agents. At lower pH values in the alkaline 
range there is complex formation with amine but as the concentration of the 
hydroxyl group becomes high as compared to the concentration of the nickel, 
hydroxyl group also enters the complex. 


4. Nature of Complex Formation in Sodium Hydroxide Solutions 


(i) Monoethanolamine.—In monoethanolamine solutions reduction 
occurs nearly reversibly. Applying Lingane’s* equation the formula for the 
complex is obtained as Ni (Moen), (OH),. 


(ii) Diethanolamine.—In this case the reduction process is irreversible. 
Applying equations for irreversible waves? one gets that one molecule of di- 
ethanolamine enters at lower concentrations of the amine and two molecules 
at higher concentrations of the amine. 


(iii) Triethanolamine.—With triethanolamine solutions the irreversibility 
of the main wave does not alter much with the change in the concentration 
of the amine in solutions containing 1 M sodium hydroxide. Application 
of equations for irreversible waves? indicates that one molecule of the tri- 
ethanolamine has entered the complex. 


5. Nature of Complex Formation in Sodium Carbonate Solutions 


(i) Monoethanolamine.—When monoethanolamine is used in the base 
solution, the reduction process is nearly reversible. From the data presented 
in Table III it is clear that one carbonate group has entered in the complex. 
Depending upon the concentration of the amine 2, 3 or 4 molecules of the 
amine enter the complex. The composition of the complexes can be written 
as follows: Ni (Moen), (COs);, Ni (Moen)3 (COs);, Ni (Moen), (CO,),. 


ible | 
of 
the 
ad 
ble 
1xi- 
ed, 
ngs 
ine 
to 
of 
na- 
to | 
ion 
led 
the 
age 
the 
um 
is | 
wo 
ili- 
ent 
the 
tial 
the 
lly 
ive 
is | 
ine 
an 
ns. | 


192 R. S. SUBRAHMANYA 

(ii) Diethanolamine.—The polarographic waves are irreversible in nature, 
It can be concluded by application of equations for irreversible waves? that 
one molecule of the amine and one CO,~ group enters the complex. The 
formula can be written as follows: Ni(Dien), (CO,),. 


(iii) Triethanolamine.—It was not possible to derive any conclusions 
regarding the nature of the complex due to progressive change of irrevers- 
bility with the complexing agent. 


6. Nature of Complex Formation in Presence of Ammonium Chloride. 
Ammonium Hydroxide 


That ethanolamines have entered into the complex is indicated by the 
shift in the half-wave potential. It can be concluded by the application of 
equations for irreversible waves” to the data presented in Table IV to the con- 


centration range where irreversibility is a constant, it is found that one mole- 
cule of the amine enters the complex. 


7. Irreversibility of Polarographic Waves 


In the case of nickel both the di- and triethanolamines give irreversible 
polarographic waves in contradistinction to the behaviour of cadmium, 
where only the triethanolamine complexes give irreversible waves. 


To explain the irreversibility of the reduction process in the case of tri- 
ethanolamine complexes of cadmium,* it has been suggested that the pre- 


sence of three hydroxyl groups in triethanolamine may be responsible to make 
the reduction process activation-controlled. 


In the case of nickel, however, the gradation of the three types of com- 


plexes is marked, the degree of irreversibility depending upon the number of 
ethanol groups in the amine. 


In the case of monoethanolamine, the reduction process is reversible. 
In base solutions containing ammonium chloride-ammonium hydroxide, 
there is a slight irreversibility. This may be due to the presence of ammonia 
in the complex. In the case of di- and triethanolamines, the reduction pro- 


cesses are irreversible, the degree of irreversibility being much greater with 
triethanolamine complexes. 


8. Height of the Polarographic Wave and the Number of Ethanol Groups in 
the Amine 


Reference to Tables II, III and IV indicates that the height of the polaro- 
graphic wave decreases with an increase in the number of ethanol groups 
indicating a decrease in the diffusion coefficient. 


9, 7 
reduc 
possi 
disso 
com| 
Con: 
reco 
cont 
Ni( 

| 
Ni 
Ni 
thi 
Te 
a 
0 
t 


Polarographic Behaviour of Metals in Ethanolamines—III 193 


9, Thermodynamic Data from Polarographic Measurements 


The reversible half-wave potential of nickel ion is not known since the 
reduction of the simple nickel ion proceeds irreversibly. It has not been 
possible to calculate the hypothetical reversible half-wave potential from the 
dissociation constant data of the ammonia complexes as the reduction of these 
complexes do not proceed reversibly at the dropping mercury electrode. 
Consequently, the conclusions arrived at will not be quite accurate. 


In conclusion, it can be stated that the polarographic investigations 
recorded in the present paper reveal the formation of the following complexes 
containing the ethanolamines: 


Ni (Moen). (OH)2 Ni (Moen), (COs); Ni (Moen); (COs), 
Ni (Moen), (COs), Ni (Moen), (NH3)z** Ni (Dien), (OH),?-¥ 
Ni (Dien), Ni (Dien), (COs); Ni (Dien), (NH3)y‘* 
Ni (Trien), (OH)z?~* Ni (Trien), (NH3)z** 


SUMMARY 


1. The polarographic behaviour of nickel has been studied in all the 
three ethanolamines at various pH values and in the following alkaline 
reagents: (1) sodium hydroxide, (2) sodium carbonate, (3) ammonium 
chloride-ammonium hydroxide. 


2. The reduction process is essentially irreversible at all pH values 
except in alkaline monoethanolamine solutions. 


3. The phenomenon of “ split-wave” is very prominent in triethanol- 
amine solutions above pH 7, and in highly alkaline diethanolamine solutions. 
An explanation based on the sluggish equilibrium between different species 
of complexes has been put forward. 


4. Evidence is presented for the existence of the following complexes 
based on polarographic data: 


Ni (Moen), (OH). Ni (Moen), (COs); Ni (Moen), (CO,); 
Ni (Moen), (CO,); Ni (Moen), (NH,)z*+ Ni (Dien), (OH)y?-¥ 
Ni (Dien), (OH),?-2 Ni (Dien), (CO,); Ni (Dien), 
Ni (Trien), (OH),?-2 Ni (Trien), (NH,),*+ 


5. Polarographic determinations of nickel can be made taking advantage 
of the reduction of the ethanolamine complexes. 
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Difficulties in calculating the thermodynamic constants from polaro. 


graphic data even in the case of reversible systems, have been pointed out, 
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EXPLANATION OF PLATES 
Fic. 1. Effect of sodium hydroxide (1-0 M) on the polarographic behaviour of nickel (1 - 600 milli- 
molar) in triethanolamine. 

Curves | to 5 have been taken using 0-02, 0-10, 0-30, 0-60 and 1-00 M triethanolamine 
respectively. There was partial precipitation under the conditions in which Curve | 
was recorded. 

Sensitivity: s/15 for Curve 1 and s/50 for Curves 2 to 5. 

The voltage line just preceding each polarogram corresponds to — 1-073 V ys. S.C. 

Fic. 2. Effect of sodium hydroxide on the polarographic behaviour of nickel (1-600 millimolar) 
in mono- and diethanolamine. 

Curves | and 2 have been obtained in base solutions containing 0-1 M sodium hydroxide 
and 0-50 and 1-00 monoethanolamine respectively. 

Curves 3 and 4 have been taken in base solutions containing 1-00 M sodium hydroxide 
and 0-50 and 1-00 M monoethanolamine respectively. 

Curves 5, 6 and 7 have been taken in base solutions containing 1-00 M sodium hydroxide 
and 0-30, 0-65 and 1-00M diethanolamine respectively. 

Sensitivity: 5/50. 

The voltage lines preceding Curves | to 4 correspond to — 0-873 V vs. S.C.E., while those 
preceding Curves 5 to 7 correspond to — 0-973 V vs. S.C.E. 

Fic. 3. Effect of sodium carbonate on the polarographic behaviour of nickel (1-600 millimolar) 


in monoethanolamine. 

Curves 1 to 3 and 4 to 6 have been taken in 0-10 and 1-00 M sodium carbonate solu- 
tions using 0-30, 0:65 and 1-00 M monoethanolamine respectively. 

Sensit vity: s/50. 

The voltage line preceding each polarogram corresponds to — 0-773 V vs. S.C.E. 
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POLAROGRAPHIC BEHAVIOUR OF METALS 
IN ETHANOLAMINES 


Part IV. Cobalt (II) 


By R. S. SUBRAHMANYA 
(Department of General Chemistry, Indian Institute of Science, Bangalore-3) 


Received February 18, 1957 
(Communicated by Prof. K. R. Krishnaswami, F.A.sc.) 


INTRODUCTION 


IN previous communications from this laboratory the polarographic behaviour 
of iron, cadmium? and nickel® have been reported. In the present investiga- 
tion a detailed study of the polarographic behaviour of cobalt (ID) has been 
undertaken to elucidate the structure of the complexes that are formed and 
find out whether the polarographic waves of the complex metal ions produced 
under different conditions are suitable for the estimation of cobalt. 


EXPERIMENTAL 


The details of apparatus and experimental procedure employed have 
been described previously.! The capillary used had a drop-time of 5 seconds 
in | N potassium chloride solution and the value of m was 1-150 mg. per 
second. The polarograms were recorded at 30° +0-1°C. 


RESULTS 


l. Effect of pH 

In highly alkaline triethanolamine solutions, no wave has been obtained. 
As the pH is lowered to 11, the reduction wave appéars but the diffusion 
current plateau is not flat due to the nearness of the fnal current rise. The 
best waves are obtained between pH 10-8 (Fig. 1). Even in these cases the 
diffusion current plateau is not sufficiently flat to enable the height of the 
wave to be accurately measured. Below pH 8, the shape of the wave is too 
poor to be of any practical use. 

In diethanolamine solutions no wave is obtained at a pH of 13-5 while 
at ca. 12-7 ill-defined wave is obtained. Lowering of pH below 11-5 improves 
the shape of the wave further. In monoethanolamine solutions, on the 
other hand, fairly well-defined polarographic waves are obtained even at 
pH 13-5, 
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The results concerning the effect of pH on the half-wave potential and 
reversibility of the wave are given in Table [. These results indicate clearly 
that the waves are mostly irreversible in nature and that the extent of 
irreversibility varies with pH. 


TABLE | 


Effect of pH on the polarographic behaviour of cobalt in presence of 
ethanolamines 


Concentration of cobalt: 1-440 millimolar; Concentration of gelatin: 0-005%; Concentra- 
tion of amine: 0-10M. 


Half-wave 
Amine used pH potential E; —E; 
E; vs. S.C.E. 
(Volt) (Volt) 
Monoethanolamine 6°5 —1-206 0-137 
9-2 —1-132 0-048 
11-0 —1-226 0-027 
12-7 —1-+340 0-040 
13-5 —1-452 0-094 
Diethanolamine .. 7:0 Ill-defined waves 
—1-209 0-101 
9-0 —1-149 0-051 
10-7 —1-216 0-034 
12-7 —1-570 0-116 
13-5 Wave not formed 
Triethanolamine .. 7-0 lll-defined waves 
7-2 —1-+293 0-075 
8-65 —1-244 0-066 
9-80 —1-279 0-044 
12-7 Wave not formed 
13-5 do. 


2. Effect of variation of the composition of the base electrolyte 


As the complex formation of cobalt takes place predominantly in the 
alkaline range, (i) sodium hydroxide, (ii) sodium carbonate, and (iii) am- 
monium chloride-ammonium hydroxide, have been employed in the base 
solutions. 


(i) Sodium hydroxide.—In Table II are presented the results of the polaro- 
graphic behaviour of ethanolamine complexes of cobalt in sodium hydroxide 
solutions. Gelatin (0-005%) was added to supress the maxima in solutions 
containing mono- and diethanolamines. 
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TABLE II 


Effect of sodium hydroxide on the polarographic behaviour of cobalt 
in ethanolamines 


Conc. of cobalt: 1-440 millimolar; Concentration of gelatin: 0-005%. 


Concentration of 


Amine used Half-wave E; —Ei 
Amine Sodium potential 
(Molar) hydroxide E; vs. S.C.E. 
(Molar) (Volt) (Volt) 
Monoethanolamine 0:10 0-10 —1-340 0-040 
0-50 0-10 —1-355 0-061 
1-00 0-10 —1-371 0-059 
0-10 1-00 —1-452 0-094 
0-50 1-00 —1-497 0-070 
1-00 1-00 —1-507 0-057 
Diethanolamine .. 0-10 0-10 —1-570 0-116 
0-50-1-00 0-10 Waves not formed 
0-10-1-00 1-00 do. 
Triethanolamine 0:10-1:00 0-10-1-00 do. 


The results of Table II show that the irreversibility of reduction of di- 
ethanolamine complex is much greater than that of the monoethanolamine. 
It is surprising that irreversibility should diminish with increasing concentra- 
tions of monoethanolamine in 1 N sodium hydroxide. 


(ii) Sodium carbonate.—Figure 2 gives some of the polarograms obtained 
with sodium carbonate solutions. It is clear from the figure that diethanol- 
amine complexes give well-defined waves while the waves obtained in tri- 
ethanolamine are complicated in nature. Curves in monoethanolamine are 
also well defined as in the case of diethanolamine. In triethanolamine there 
are two or three waves very close to each other. Even in the case of diethanol- 
amine there is a slightly elongated bend before the diffusion current plateau. 
The results of analysis of the waves are given in Table III. 


(iii) Ammonium chloride and ammonium hydroxide.—Figure 3 shows 
the polarograms obtained in presence of diethanolamine. Curves obtained 
in mono- and triethanolamines are similar in nature. Table IV indicates 
the results of the analysis of the waves. From the table it is clear that the 
waves are slightly irreversible in all the three amines. 
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TABLE [Il 


Effect of sodium carbonate on the polarographic behaviour of cobalt 
in the presence of ethanolamines 
Concentration of cobalt: 1-440 millimolar; Concentration of gelatin: 0-005%. 


; Concentration of 
Amine used Half-wave 


Amine Sodium potential E; —E 
carbonate BE; vs. S.C.E. 
(Molar) (Molar) (Volt) (Volt) 
Monoethanolamine 0-10 0-10 —1-256 0- 018 
0-50 0-10 —1-273 0-028 
1-00 0-10 —1-289 0-026 
0-10 1-00 —1-287 0-027 
0-50 1-00 —1-297 0-032 
1-00 1-00 —1-314 0-026 
Diethanolamine 0-10 —1-298 0-066 
0-50 0-10 —1-374 0-109 
—1-58 
1-00 0-10 —1-429 0-105 
—1-53 
0-10 1-00 —1-333 0-070 
0-50 1-00 —1-437 0-122 
1-00 1-00 —1-508 0-128 
Triethanolamine .. Two or three poorly separated and ill-defined polagrams 


not amenable to analysis. 


TABLE [V 


Effect of ammonium chloride (1 M)—ammonium hydroxide (\ M) 
on the polarographic behaviour of cobalt in ethanolamines 
0-005% geatin used as maximum suppressor. 


Concentration Half-wave 
Amine used of amine potential E; —E; 
Ey vs. S.C.E. 
(Molar) (Volt) 
0-10 —1- 2712 0-044 
0:50 —1-293 0-042 
1-00 —1-305 0-038 
Diethanolamine 0-10 —1-279 0-053 
0:50 —1-300 0-048 
1:00 —1-311 0-040 
Triethanolamine 0-10 —1-281 0-044 
0-50 —1-+321 0-049 
1-00 —1-339 0-050 
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3. Effect of concentration of cobalt on the diffusion current constant in 
presence of sodium hydroxide, sodium carbonate and ammonium chloride 
(1 M)—ammonium hydroxide (1 M) 


The results presented in Tables II and III indicate that the polarographic 
waves in sodium hydroxide solutions containing diethanolamine and tri- 
ethanolamine are complicated and ill-defined and hence unsuited for the 
polarographic estimation of the metal. For a similar reason solutions con- 
taining sodium carbonate and triethanolamine are also not suitable. The 
effect of concentration of cobalt (II) on the diffusion current constant has 
been tested in base solutions containing monoethanolamine + sodium hydro- 
xide, monoethanolamine + sodium carbonate, and diethanolamine + sodium 
carbonate where we get very well-defined waves. 0-005% or 0:025% gelatin 
has been included in the base solutions to eliminate the maximum. Due to the 
limited solubility of the complexes at low concentration of the amine, 1-00 M 
monoethanolamine and 0-5 M diethanolamine solutions have been used. Higher 
concentration of diethanolamine is not of much advantage as the diffusion 
current plateau is not well defined. As the addition of ethanolamines in 
base solution containing ammonium chloride-ammonium hydroxide does not 
have any practical advantage in the estimation of cobalt (ID, these solutions 
have been omitted in testing the effect of concentration of cobalt (II) on the 
diffusion current constant. The results obtained are given in Table V, from 
which it is clear that these base solutions are suitable for the estimation of 
cobalt (I). 


DIsCussION 
(i) Effect of pH 
The results given in Table I indicate that when the pH of the solution 
is decreased from 9 to about 7, the half-wave potential and the value of E;—E, 
increase. The half-wave potential at pH 7 is — 1-21 V in mono- and di- 
ethanolamine and — 1-29 V in the case of triethanolamine. The above 
results indicate that the reduction step may be due to 


Co (H;,O),** — Co. 
This, however, does not seem consistent with the fact that the half-wave 
potential of hexacobaltous ion is — 1-425 V us. S.C.E. A similar anomaly 
has also been noticed by Willis, Friend and Mellor,* who noticed a Ej of 


— 1-28 V for the second wave of the Co**+ complexes. This wave has been 
interpreted by the authors as due to the reduction of Co(H,O),*+ + Co. 


When the pH of the solution is increased to 9, the half-wave potential 
decreases and reduction process becomes more reversible in character. This 
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TABLE V 


Effect of concentration of cobalt on diffusion current constant* 


Concentration Diffusion Concentration Diffusion 
of cobalt current of cobalt current 
(Millimolar) constant (Millimolar) constant 
ig ig 
Cmi tt Cmi tt 
(1) B.S.: 1-00 M Moen-+-0-10 M (4) B.S.: 1:00 M Moen+1-00 M 
NaOH+0:-005% gelatin Na,CO,+0-025% gelatin 
X = 1-381 at — 1-50 V vs. S.C.E. X=1-381 at —1-60 V vs. S.C.E. 
1-440 3-04 1-440 2:55 
7-200 3-05 7-200 2°57 
14-400 3-05 14-40 2-48 
(5) B.S.: 0°50 M Dien+0-10 M 
(2) B.S.: 1-00 M Moen+1-00 M Na,CO,+0-005% gelatin 
NaOH-+-0-005% gelatin X=1-360 at —1-75 V. vs. S.C.E. 
1-440 2-86 1-440 2-98 
7-200 2:85 7-200 2°95 
14-400 2-84 14-40 2:87 
(6) B.S.: 0:50 M Dien+1-0 M 
(3) B.S.: 1-00 M Moen+0-10 M Na,CO,+0-005% gelatin 
Na,CO,+0:-025% gelatin X=1:-345 at —1-75 V vs. S.C.E. 
X=1-381 at —1-60 V. vs. S.C.E. 1-440 2-38 
1-440 3-12 7-200 2-36 
7-200 3-06 14-40 2-38 


* Abbreviations—B.S.—Base solution; X—m?/*f/6 in mg?/? sec..-!/2;_ Moen—monoethanol- 
amine; Dien—diethanolamine: Trien—triethanolamine. 
indicates complex formation with the amine. At still higher pH values the 
hydroxyl group enters into the complex forming hydroxo-ethanolamino 


complex. The second increase in E; above pH 11 is perhaps due to the entry 
of OH- group into the complex. 


(ii) Splitting of the polarographic wave 


Sodium carbonate solutions—The reduction of Co(II) at the surface 
of the dropping mercury electrode is a stright two-electron step. Hence the 
splitting of the polarographic wave observed in triethanolamine solutions 
is due to two or three different species? in sluggish equilibrium. Similar 
explanation can be given for the elongated bend that appears before the 
diffusion current plateau in diethanolamine solutions. In monoethanolamine 
solutions, the splitting is practically absent. It is noteworthy that the ten- 


dency for splitting decreases progressively with the decrease in the number 
of ethanol groups in the amine. 
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Ammonium chloride-Ammonium hydroxide solutions——The slight ten- 
dency for the splitting of the wave as indicated by a pressed-in bend before 
the diffusion current plateau is evident only when triethanolamine is used 
as the complexing agent. 


(iii) Irreversibility of the polarographic waves 


A perusal of the Tables I, II, III and IV indicates that the reduction of 
the ethanolamine complexes is essentially an irreversible process, excepting 
in the case of the monoethanolamine complex in presence of sodium carbonate 
solutions. In contrast to the behaviour of nickel* complexes, even some of 
the monoethanolamine complexes of cobalt are also irreversible and the 
extent of irreversibility increases with the number of ethanol groups in the 
amine. It is surprising that the hydroxomonoethanalo complex should be 
reduced irreversibly; perhaps, the reduction of the complex is activation 
controlled under these conditions. 


(iv) Nature of complex formation in sodium hydroxide solutions 


Application of the modified equations for irreversible processes! to the 
results presented in Table II indicates that one molecule of monoethanolamine 
is present in the complex. It is not possible to compute the number of 
hydroxyl groups in the complex as the reduction processes in 0-1 Mand 1:0M 
sodium hydroxide differ in irreversibility. 


(v) Nature of complex formation in sodium carbonate solutions 
Application of Lingane’s equation for reversible processes? indicates 
the following formula for the complex: 
Co (Moen), (COs). 


For reasons similar to those explained in the previous paper,* it has not been 
possible to calculate the thermodynamic constants of this complex. 


In the case of diethanolamine, the reduction process is irreversible. From 
the equations for irreversible processes! it follows that two diethanolamine 
groups and one carbonate group enter the complex: 


Co (Dien), (CO3);. 
(vi) Nature of complex formation in presence of ammonium chloride and am- 
monium hydroxide 


From the results presented in Table IV it is clear that the ethanolamines 
enter into the ammonia complex. Applying the equations for irreversible 
processes,! it is found that one molecule of ethanolamine has entered into 
the complex. 
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In conclusion it may be stated that the polarographic studies of the 
cobalt complexes of ethanolamines have revealed the formation of the fol- 
lowing complexes : 


(1) [Co (Moen), (OH),,]?-* (2) Co (Moen), (COs); 
(3) Co (Dien), (COs), (4) Co (Moen), (NH3)y** 
(5) Co (Dien), (NH3)2** (6) Co (Trien), (NH3)q**. 
It may also be pointed out that the reduction waves obtained in mono- 
and diethanolamines can be used for the estimation of cobalt. 
SUMMARY 


1. The polarographic behaviour of cobalt has been studied in ethanol- 
amines employing the following electrolytes to alter the pH of the solutions: 
(i) Sodium acetate-acetic acid, (ii) Sodium carbonate, (iii) Sodium hydroxide, 
and (iv) Ammonium chloride-Ammonium hydroxide. 


2. Studies of the effect of pH have indicated that in acidic solutions 
aquo-cobalt ion is reduced at the dropping mercury electrode. As the pH 
is increased to 9, the amine complex is formed, and when the pH is increased 
still further, hydroxyl group also enters into the complex. 


3. The reduction of the ethanolamine complexes of cobalt at the dropping 
mercury electrode is essentially an irreversible process except in solutions 
containing sodium carbonate and monoethanolamine. 


4. The foregoing studies indicate the formation of the following com- 
plexes of cobalt :— 
(1) [Co (Moen), (OH)z]?~* (2) Co (Moen), (COs); 
(3) Co (Dien), (COs), (4) Co (Moen), (NHs3)y** 
(5) Co (Dien), (NH;)2** (6) Co (Trien), (NH3)q**. 
5. It has been indicated that the reduction waves of cobalt in mono- 
and diethanolamines can be used for the analytical estimation of the metal. 
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EXPLANATION OF PLATES 


Effect of pH on the polarographic behaviour of cobalt (1-440 millimolar) in triethanol- 
amine. 


Curves 1 and 2 have been taken at pH 9-8 and 8-65 respectively. 


Curves 3 and 4 are for the reduction of Co+++, — Cot+, Cot++ being formed by the 
oxygen present in the solution. 


Sensitivity: 5/30. 

The voltage lines just preceding Curves 1 and 2 correspond to — 1-073 V vs. S.C.E. 
and those preceding Curves 3 and 4 correspond to 0-027 V vs. S.C.E. 

Effect of sodium carbonate (0-1 M) on the polarographic behaviour of cobalt (1-440 milli- 
molar) and di- and triethanolamine. 

Curves | to 3 are obtained in 0-10, 0-50 and 1-00 M diethanolamine respectively. 

Curves 4 to 7 are obtained in 0-10, 0-50, 1-00 and 0-02 M triethanolamine respectively. 

Sensitivity: 5/50. 

The voltage lines just preceding polarograms | to 3 correspond to — 1-073 V vs. S.C.E. 
and those preceding Curves 4 to 7 correspond to 1-273 V vs. S.C.E. 


Effect of ammonium chloride (1-0 M)—ammonium hydroxide (1-0 M) on the polaro- 

graphic behaviour of cobalt (1-440 millimolar) in presence of diethanolamine. 
Curves | to 3 are obtained employing 0-10, 0-50 and 1-00 M diethanolamine respectively. 
Sensitivity: 5/50. 


The voltage lines just preceding the polarograms correspond to — 1-073, —1-:073 and 
— 1-173 V vs. S.C.E. 
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ERRATA 


Vol. XLV, No. 2, Section A, February 1957 
Line 
or For 
Equation 
6 ou ou 
ror roo 


333 —K, +Ky 


4.1, 4.6 +, 


5.4, 5.5, 5.6, 5.7 —K, +Ky 
5.9 Delete — 2 Ky (1 — c) 


5.13 and last equation —K,y +Kgy 
5.11 Multiply the right-hand side, giving 7 
the values of 

a’ and C,? by (3 — 2c) 
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NOTICE TO AUTHORS 


Scientific papers intended for publication in the Proceedings of the Indian 
Academy of Sciences can be accepted only when they are communicated by a 
Fellow of the Academy whose duty shall be to satisfy himself that such communi- 
cations are fit to be read at the Meeting of the Academy and published in: its 
Proceedings. 

Papers should not ordinarily exceed fifty pages of foolscap. MSS. should 
be either typewritten or written in legible hand on one side of the paper. All 
papers should be carefully revised by the authors and should be absolutely in final 
form for printing. Position for text-figures should be indicated. Each paper 
shall conclude with a critical summary not exceeding 350 words. 

Drawings, diagrams or other illustrations should be made on larger scale 
(preferably) twice the size than the ones in which they are intended to appear. 
They should be done in Indian ink on bristol board with lettering in pencil. Scale 
of magnification of camera lucida tracings should be indicated by the side of 
drawings. In certain special cases arrangements will also be made for mono- 
chrome lithographic and other colour plates. Reduction of illustrations desired 
should be indicated in pencil. Appropriate legends should accompany all drawings. 
Names of authors are to be marked in pencil on the left-hand corner of drawing 
sheets. Photomicrographs should be securely mounted with colourless paste. 

All tables, quotations and footnotes which will be set hereafter (beginning 
from Vol. I, No. 2) in types smaller than the text, should be typewritten on sepa- 
rate sheets and placed with the text in proper sequence. Footnotes should be 
numbered in Arabic numerals. 

References to literature in the text should be given, whenever possible, in 
chronological order, only the names of authors and years of publication, in 
brackets, being given. They should be cited in full after the summary, the 
authors’ names following in alphabetical order. Thus, 

Name or Names of author; Name of Journal (abbreviation) with a single 
underline; Year of publication; Number of Volume with a double underline, and 
lastly page. The following would be a useful illustration:— 

Bergmann and Stather Z. Physiol. Chem., 1926, 152, 189. 

Two copies of slip-proof and wherever possible, a page proof for final 
revision will be sent to authors. All corrections are best made on the slip-proof 
which should be transmitted to the Office of the Academy. All proof corrections 
involve heavy expenses which would be negligible if the papers are carefully 
revised by the authors before submission. 


Fifty free reprints including plates and with cover will be supplied for each 
paper. Additional copies can be supplied at cost on previous intimation. 

Blocks appearing in the Proceedings will be available for purchase by their 
respective authors. Orders for the same should be sent along with the corrected 
proofs and in any case not later than one month after the date of publication 
of the paper. The price charged would be 25% of the actual cost of the blocks 
plus freight and despatching charges. If the blocks are reproduced in other 


journals or publications, due acknowledgment should be made in them to the 
Proceedings. 


The original drawings and plates of blocks appearing in the Proceedings 
will be returned to such of the authors as may require them provided the cost of 
despatching such originals is borne by them. 
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